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ABSTRACT 

We present direct imaging observations at wavelengths of 3.3, 3.8 (L'band), and 4.8 (M band) /im, 
for the planetary system surrounding HR 8799. All three planets are detected at L' . The c and d 
component are detected at 3.3 ^tm, and upper limits are derived from the M band observations. These 
observations provide useful constraints on warm giant planet atmospheres. We discuss the current age 
constraints on the HR 8799 system, and show that several potential co-eval objects can be excluded 
from being co-moving with the star. Comparison of the photometry is made to models for giant 
planet atmospheres. Models which include non-equilibrium chemistry provide a reasonable match to 
the colors of c and d. From the observed colors in the thermal infrared we estimate T e ff< 960 K for 
b, and T e ff=1300 and 1170 K for c and d, respectively. This provides an indepen dent check on the 
effective temperatures and thus masses of the objects from the iMarois et al.l (|2008|) results. 
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1. INTRODUCTION 

Over the past decade, a number of techniques have 
dramatically expanded our understanding of the nature 
of exoplanets. Initial detection of systems via radial 
velocity variations has been followed up by studies of 
transits, astrometric confirmation, and gravitational mi- 
crolensing. All of these approaches have been helpful in 
developing our picture of planetary system architectures, 
and providing insight into formation. 

The direct imaging of extrasolar planets is the latest 
technique to provide useful information, with detection, 
of ob jects, first, around low-mass objects dChauvin et al.l 
2004), followed, more recently, by the dete ction of sev- 
eral planets around intermediate mass stars (|Kalas et al.l 
120081: IMarois et al.ll2008t iLagrange et al.ll2008fl ~ 

Direct images of extrasolar planets not dominated by 
insolation have the potential to provide a wealth of in- 
formation about the size, temperature, composition, and 
even formation history of these objects. Th e recently dis- 
cover ed planetary system around HR 8799 (jMarois et al.1 
2008) , with three massive planets at large orbital separa- 
tions, may provide one of the most useful laboratories to 
constrain the spectral energy distribution of young giant 
planets. 

HR 8799 is a young A5V star, thought to b e ap- 
proximately 30-160 Myr old (jMarois et all 120061 ). It 
is kn own to have a bright debr is disk at large orbital 
radii (' Williams fe Andre ws 2006) with approximately 0.1 



thought to exist as well (|Chen et all l2006h . Taken to- 
gether with the planets, this suggests a system with an 
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inner dust disk truncated by the d component at 24 AU, 
and an outer disk truncated by the b component beyond 
~ 80 AU. Su et al. (2009) also find an outflow of small 
grains, likely caused by gravitational stirring from the 
planets. The star also has a number of other interesting 
properties, including X-ra y emission (|Heartv et allll999t 
iSchroder fe Schmittll2007T ) and a deficiency in refractory 
metals in the stellar atmosphere (the A Boo phenomenon; 
iGrav fe Kayll99ftlGray et alJl2003f ). Multiple metallic- 
ity esti mates are consisten t it being metal-poor (Fe/H 
~ -0.5: lGray fc Kayelll999( l but having solar abundance 
lighter elements, such as c arbon and oxygen ([Sadakand 
12005 iGerbaldi et al.ll2007l ). 

The planets around HR 8799 are interesting, particu- 
larly in the context of planet formation alternatives. If 
the planets formed in s itu, the core accretion hypoth- 
esis ([Pollack et al.l Il996h would require formation of 10 
M^orfd cores at distances of 40 and 70 AU well befo re the 
dispe r sal of the gas. A disk fragmentation scenario dBoss 
19971; iNero fc Biorkmanl l2009j : iDodson-Robinson et"al 



20091 iHelled fc Bodenheimerl l2009) may be able to more 
naturally explain the massive planets and their location. 
A more detailed look at the planet's environment and 
their spectral energy distributions (SEDs) can provide 
clues to these alternatives. 

M odels of the spectral energy distribution for plan- 
ets (|Burrows et al.l 119971: iBaraffe et all 120031 ) are used 
to estimate the temperature and mass of the objects. 
However, these models are currently only constrained 
by field brown dwarfs and objects such as 2M1207 , 
whose formation history is uncertain (|Lodato et alj|20"0"5t 
lMamaiek""fe Meyer 2007]). Are giant planets formed in a 
circumstellar disk around a normal star different? Ob- 
jects such as the planets ar ound HR 8799, as well as the 
planets around Fomal haut (jKalas et al.ll2008| ) and j3 Pic 
(jLagrange et al.ll2008l ) may provide our first opportunity 
to address this topic. With multi-wavelength constraints 
of the brightness of the planets we can disentangle effects 
from clouds, composition, and vertical mixing that may 
affect the measurement of the temperature, and thus the 
mass of the planets. By comparison of observations with 
field objects and current models, we may be able to de- 
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velop a better understanding of the physics leading to 
the SED's of giant exoplanets. 

The 3-5 (an region provides access to both the 3.4 
/im CH4 feature and the CO bandhead at 4.7 /im. Ob- 
servations of this portion of the SEDs provides useful 
constraints on the relative amount of CO to CH4 in the 
atmospheres of these planets. The relative absorption of 
these species, in the near-infrared, is used to define the 
transition in spectral type between the L and T sequence 
among field brown dwarfs. 

Observations of brown dwarfs in this region 
(jLeggett et al.l 120071 ) indicate that their colors are 
best reproduced by models that have substantial vertical 
mixing between the hot lower layers and the cooler 
upper atmosphere. Comparison of HR 8799 planets to 
these objects will allow us to better understand how 
closely exoplanet atmospheres may mimic those of field 
brown dwarfs. 

Constraining the SED of giant planets in this region 
may help guid e future direct ima ging planet searches. 
Both Jupiter dGillett et all I1969D and Gliese 229 B 
(jOppenheimer et al.lll995l ) have a broad peak in the flux 
at 4-5 /im, suggesting that this may be a robust fea- 
ture of brown dwarfs and giant planets in the T e ff=100- 
f 000 K range. However, variou s models indicate diff erent 
bands as being prefer able (c.f. iMarlev et al.l (|1996l ). and 
iBurrows et al.l (|1997l )). In order to understand whether 
searching in L' or M band is preferable, it would be help- 
ful to understand giant planet colors over the expected 
temperature range of these objects. 

In this paper we present observations at 3.3, 3.8 and 
4.8 /im of the planets orbiting HR 8799. In Section 2 we 
describe the MMTAO and Clio camera observations of 
the HR 8799 system. Section 3 details the data reduc- 
tion and analysis. Section 4 describes the photometric 
and astrometric results. Section 5 analyzes the age es- 
timates for HR 8799. Section 6 compares these results 
to theoretical models and observations of brown dwarfs 
in the field. Finally, in Sections 7 and 8 we discuss the 
interpretation of the results and summarize our conclu- 
sions. 

2. OBSERVATIONS 

HR 8799 was observed on 21 November 2008, 8 
January 2009, and 12 September 2009, using the 
MMT deformable secondary-based adap tive optics sys- 
tem (jWildi et al.M2003t iBrusa et al.ll2004f ). Cirrus clouds 
intermittently affected the throughput of the observa- 
tions on 21 Novem ber 2008, and 12 September 2009. 

Th e Clio camera (|Freed et al.ll2004J : iSivanandam et all 
I2006D. an optimi z ed instrument for exoplanet detection 
(jHinz et al.ll2006l : iHeinze et al.l l2006) in this spectral re- 
gion, was used for the observations. The camera has a 
high flux detector that allows high duty cycle observa- 
tions in the L' and M band atmospheric windows. Com- 
bined with the low backgroun d provided by a deformable 
secondary (|Llovd-Hartl |2006T) . the system is capable of 
sensitive observations in the thermal infrared. 

The Clio observ ations use the Barr M ban d (A=4.8 
jum, AA=0.6 /im, ISimons fc Tokunagal {2002)), rather 
than the more common MKO M' passband (A=4.68 /im, 
AA=0.22 /im). The Barr M band filter is preferred 
for sensitive M band observations that maximize the 
signal-to noise for the expected SED of a cool object 



(|Freed et al.ll200l . At L' the MKO filter was used. 
The zeropoints used for these filters are 248 Jy and 154 
Jy for L' and M resp ectively (|Tokunaga &; Vaccal 120051 : 
IBessell fc Br^[l988ft . 

In addition to standard MKO L' and Barr M band fil- 
ters, we carried out observations with a filter which was 
centered at 3.3/xm with half-power points at 3.10 /im 
and 3.5 /im. This filter is well matched to the methane 
absorption located at 3.4 /im and lies in the short wave- 
length region of the L band atmospheric window. How- 
ever, it is not closely related to any standard photometric 
system in the L window. We refer to it in the rest of the 
paper as the 3.3 /im filter, and list magnitudes for this 
filter with the notation [3.3]. The zeropoint adopted for 
this filter, is 330 Jy, obtained by interpolating between 
the L' and K band photometric zeropoint. 

Clio observations are obtained by acquiring individual 
integrations of the detector with frame times set to avoid 
saturation by the background emission, typically 0.2 s for 
M band, 1.5 s for L' band. These individual exposures 
are accumulated into cubes of data with approximately 
15-20 seconds total integration. Typically, five exposures 
are obtained before nodding the telescope in elevation. 
This results in the star being moved approximately ev- 
ery minute from one side of the detector to the other to 
provide periodic measurements of the sky background, 
as well as detector variations. Each nod pair provides a 
complete set of data with ten independent measurements 
of the field around HR 8799. 

A significant source of noise in the images are static or 
slowly changing speckles. The effect of these are mini- 
mize d with the technique of angular d ifferential imaging 
(e.g. lMarois et aLlteOOalHinz et al.l200H . During the ob- 
servations the instrument derotator is not operated. This 
has the effect of causing the sky coordinates to slowly ro- 
tate relative to the detector coordinates. Since the static 
speckles stay fixed with respect to the telescope, we can 
remove them in the final image processing without af- 
fecting the flux from an off-axis companion. 

Observations were taken in both L' and M band on 
November 21, 2008. Approximately 70 degrees of rota- 
tion was obtained at L' and 30 degrees at M band (see 
Table [T]). For L' band observations, 30 frames were re- 
jected due to image jitter while observing near transit. 
The final image (see Figure [l| uses approximately 90 
minutes of data. The M band observations were taken 
just after this, as intermittent cirrus clouds began to 
form. Consequently, although approximately 170 min- 
utes of integration was acquired, only 84 minutes of total 
integration (210 total frames) were used, after rejecting 
those images which showed stellar fluxes less than ap- 
proximately 60% of the "clear" images. 

A first set of 3.3 /im observations were taken on Jan- 
uary 9 2009, which required observing the star well af- 
ter it had transited the meridian. This resulted in only 
a small amount of sky rotation for these observations. 
Consequently, for the 3.3 /im observations we observed a 
standard star for use in PSF subtraction. A total of 48 
minutes of integration was obtained in the 3.3 /im filter. 

A second set of 3.3 /im observations were taken on 12 
September 2009. A total of 113 minutes of integration 
was obtained, with approximately 30 degrees of rota- 
tion during the observations. Light clouds were present 
throughout the observations. 
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As can be seen in Figure [TJ slowly changing speckles 
have different levels of importance for the three observed 
wavelengths. At 3.3 and 3.8 /xm the noise in the resid- 
ual subtraction begin to degrade the detcctability inside 
approximately an arcsecond. However, for 4.8 /xm ob- 
servations, the brighter sky background and more stable 
PSF at longer wavelengths result in the images being 
sky-background limited even for the closest planet sepa- 
ration. 

3. DATA REDUCTION AND ANALYSIS 

The set of data was reduced with a custom-developed 
pipeline using the Perl Data Language (PDlQ). A cube 
of frames for each observation is coadded to obtain in- 
dividual images of 15.3, 14.6, or 24.0 s for the 3.3 /xm, 
L' and M band observations respectively. Each image 
is subtracted from its pair in the other nod position to 
create two sky-subtracted images of the star. To elimi- 
nate the effect of slightly different gains among the four 
different channels, a blank portion of the detector is used 
to measure the average flux in each channel. This value 
is subtracted from all pixels corresponding to the appro- 
priate channel. Finally, a bad pixel rejection routine is 
carried out which looks at the detector response for two 
successive nods to identify pixels which vary by more 
than 3 times the standard deviation. The values for these 
pixels are replaced by an average of their neighbors. The 
images are then shifted by measuring the centroid of the 
star and rotated by the Parallactic angle of the obser- 
vations to create a final image of the star with North 
up and East to the left in the image. A median PSF 
image is created by median combining all of the images 
without rotation. By subtracting this median PSF from 
all of the images prior to rotation, the majority of the 
diffraction and static aberration-induced scattered light 
can be removed from the final, combined image. 

The measurements of the brightnesses of the planets 
are calculated relative to the measured flux from HR 
8799 itself. Since the star is saturated in the L' and 
3.3 zxm images, images taken with neutral density filters 
or shorter exposures were obtained to measure the flux. 
The photometry is scaled from the 2MASS K s band flux 
of HR 8799, using the ex pected K-L' and K-M colors fo r 
HR 8799's spectral type (jMarois et alj|200l [Co^[2000f) . 
see Table El 

The scale of the Clio pixels has been measured to be 
0.04857 ± 0.0000 3 "/pixels, u sing periodic observations 
of binary stars (|Heinzd I2007D . The orientation of the 
detector relative to the sky is also calibrated with these 
measurements. North on the detector is 92.27±0.2 de- 
grees counter clockwise from the top of the array. 

4. OBSERVATIONAL RESULTS 

To calculate the flux of the planets, apertures of sizes 
of 0.25 arcsec were placed at the approximate locations of 
the sources. The images of the stars were measured with 
the same aperture to account for aperture-dependent ef- 
fects. To estimate the errors in these measurements, each 
data set was split in half. For each half of the data set 
a measurement of the image in each nod position was 
carried out, resulting in four independent measurements. 
The apparent magnitude and standard error in the mean 

7 http://pdl.perl.org 



for each set are reported in Table [31 For M band, and 
some of the 3.3 /xm data the planets are not detected. 
For these data, 3 a upper limits are derived by measur- 
ing variations in flux with the same size aperture, at the 
same radial separation as the planet in an azimuthal arc 
around the star. We note that, typically, we conserva- 
tively quote 5-7 a limits for detection with Clio observa- 
tions. However, for these observations, where we know 
the exact location of the planet, 3 a limits are more real- 
istic estimates of our uncertainty in the flux contribution. 

For comparison with models and other p hotometry we 
adop t a distance of 39.4 pc to HR 8799 (jvan Leeuwenl 
2007;. The photometric measurements of planets b, c, 
and d are reported in Table [3] We tabulate absolute and 
apparent magnitudes, as well as the "absolute" spectral 
flux density in milliJanskys, or the flux at a distance of 
10 pc. This is useful for comparisons with model spec- 
tral energy distrib utions. The L' ban d measurements 
are consistent with iMarois et al.l (|2008l ) , within our pho- 
tometric uncertainty. 

The initial observations at [3.3] detected c, but did not 
detect the other two components. Follow-up observations 
in September 2009 confirmed the detection of c, and also 
detected d. We report the photometry for 3.3 /imfrom 
the September 2009, as it is more sensitive than the Jan- 
uary 2009 data, with better PSF subtraction. For the 
September 2009 data, there is a marginal source at the 
location of b. However, it is below the significance level 
of the image. We thus adopt the limiting magnitude for 
the image as the upper limit for b at 3.3 /im. 

The astrometric measurements are calculated with a 
" center-of-light" algorithm within the photometric aper- 
ture. The same procedure of dividing up the data into 
four independent sets is used to calculate the astromet- 
ric uncertainty. These uncertainties are consistent with 
the astrometric measurement errors expected due to the 
PSF width of these observa tions and their low signal-to- 
noise ratio (|Lindegrenl fr978') . Astrometric results for the 
observations are summarized in Table [4] The L' position 
of planets b and c are consistent with the Marois results, 
differing by 24 and 44 mas respectively. The L' d po- 
sition is displaced by 62 mas from the reported Marois 
position. Since the orbital motion expected for this ob- 
ject is only 7 mas over the two month time span between 
the two observations, the discrepancy of d is not easy to 
explain. It is possible that the positional difference is 
an indication of contamination of the flux by a residual 
speckle in the final image. The [3.3] position of d is also 
66 mas from the Marois position, but since the observa- 
tions were taken a year after those results, the difference 
is similar to the expected orbital motion. This effect can 
be seen in Figure [TJ where the position of d in the [3.3] 
image is significantly shifted relative the the L' image. 

5. ON THE AGE AND STELLAR MULTIPLICITY OF HR 

8799 

Understanding the age of the HR 8799 system is 
important for interpreting the SEDs of the planetary 
companions and assessing their sta bility. For example 
iFabrvckv fc Murrav-Clavi (|2008l ) and lReidemeister et all 
(|2009f ) analyze the stability of this s ystem, for wh i ch one 
of the crucial parameters its age. iMarois et al.l (2008) 
assign an age for HR 8 799 o f 60±|g° (30-160) Myr, and 
previously iMoor et al.1 (2006) estimated an age of 20-150 
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Myr based on ass igning HR 87 9 9 mem bership to the Lo- 
cal Association. iMarois et al.l (|2008l ) based the age on 
careful consideration of four diagnostics: 1) the star's 
space motion and possible association with HD 984 and 
HD 21318, as well as the Columba and Carina Associ- 
ations, 2) its position in a color-magnitude diagram, 3) 
the typical age of A Boo and 7 Dor variables, and 4) the 
large mass of its debris disk. Of these, only the star's 
position on a color-magnitude diagram provides a quan- 
titative constraint. Here, we explore whether association 
of HR 8799 with any additional companions or moving 
groups can be used to further refine its age estimate. 

5.1. On the Association with HD 984 and HD 21318 

IMarois et al.l ()2008D calculate a space motion for HR 
8799 of UVW = (-11.9, -21.0, -6.8) km/s, and suggest 
that it is similar to other young stellar groups in the so- 
lar neighborhood (e.g. Columba, Carina, etc.) as well as 
two neighboring stars: HD 984 and HD 21318 to which 
they assign ages of ~30 and ~100 Myr, respectively. We 
investigate the space motion of HR8799 here and discuss 
Columba and Carina groups in the followi ng subsection 
Using the revised Hipparcos astrometry Ivan Leeuwenl 
( 2007) and the radia l veloc ity from the compiled RV cat- 
alog of iGontcharovl (|2006f) . we estimate the velocity of 
HR 8799 to be UVW = (-12.2, -21 2, -7.2) ± (0 . 6, 0.8 , 
1.2) km/s. Using astrom etry from Ivan Leeuwenl (|2007ft 
and radial velocities from lHolmberg et al.l (|2007F . we in- 
dependently verify that the stars HD 984 and HD 21318 
have space motions within 3.6 ±2.0 km/s and 1.0 ±1.8 
km/s of HR 8799, respectively. The HD 984 and HD 
21318 are separated by 26 and 19 pc from HR 8799, re- 
spectively. We simulated the orbits of HR 8799 and these 
two stars using the ep icyclic approximation, ado pting the 
Oort constants fr om iFeast fc Whitelockl ()1997D and the 
LSR veloc ity fromlDehnenl (119981) . and the local disk den- 
sity from Ivan Leeuwenl (|2007ft . We find that HD 8799 
and HD 984 were not significantly closer to one another 
within the past 100 Myr (although some ~0.02 pc closer 
~0.5 Myr ago), and were probably separated by ~80 pc 
~30 Myr ago. We find that HR 8799 and HD 21318 
would have been separated by ^140 pc ^100 Myr ago, 
with the nearest pass being ~13 pc away ~11 Myr ago. 
While ~1 km/s uncertainties in the velocities can lead to 
large uncertainties in the distant past (~1 pc/Myr), any 
genetic tie between HR 8799 and HD 984 and HD 21318 
would seem tenuous at best. 

5.2. On the Association with Columba and Carina 

Is HR 8799 related to the putative Columba or Carina 
groups? We adopt our revised space motion for HR 8799 
and the parameters of the veloc ity and posit i on cen troid 
for Columba and Carina from iTorres et al.1 (|2008f ) and 
ran the stars through our epicyclic orbit code. HR 8799 is 
cu rrently 98 pc f rom th e centroid of Columba, as defined 
by ITorres et al.l (|2008f) , and its closest approach to the 
group within the past 160 Myr was a pass ~27 Myr ago 
(A ~ 58 pc), which seems too distant to be a plausible 
candidate for ejection/evaporation. Similarly, HR 8799 is 
currently ~116 pc distant from the centroid for Carina, 
and its closest approach within the past 160 Myr was 
~37 Myr ago (A ~ 45 pc) - also relatively distant for 
a plausible association. We find that the centroids of 



Columba and Carina were closest to one another ~18 
Myr ago with a separation of ^24 pc. While the putative 
groups themselves are very large (tens to >100 pc), it is 
more likely that they constitute "complexes" with a wide 
range of ages, making membership - even if likely - not 
useful for assigning ages. 

5.3. On the Candidate Common Proper Motion Stars 

Found by Close & Males (2009) 

IClose fclV lalcs (2009) found no evidence for bound 
companions using archival Gemini and HST/NICMOS 
images. They also identified a few candidate com- 
mon proper motion stars withi n a degree using 
NOMAD /USNO-B1.0 astrometry (IMonet et al.1 l200a 
iZacharias et al.l 12005!) . however none are particular 
strong candidates for being associated with HR 8799. 
Despite the similarity of their proper motions to that 
of HR 8799, the two top co-moving c ompanion candi- 
dates identified bv lClose fc Males! ((2009) can be ruled out 
through further consideration of their photometry. NO- 
MAD #1115-0634383 (J2000; 23 10 22.93 +21 34 20.7) 
is 48' distant from HR 8799 and is very red. We find 
that the 2MASS JHK photometry for this object (J-H 
= 0.58 ±0.04 mag, H-K s = 0.22 ±0.04 mag) is consis- 
tent with an unreddened M2V (±1 subt ype). This cor- 
respon ds well with the B-V estimated bv lClose fc Males! 
(2009) (B-V ~ 1.61; inferred from NOMAD photographic 
magnitudes) and with a V-K s color (~ 4.5 mag) that we 
calculate using USNO-A2.0 B & R photometry and 
2MASS K s . However, if the star is an unreddened M1V- 
M3V dwarf at d = 39 pc (My ~ 13.9 mag), it is ap- 
proximately 2-4 magnitudes below the main sequence, 
and hence could not be related to HR 8799. Similarly, 
the second best candidate companion NOMAD #1108- 
0634609 (TYC 1717-1120-1) can be r ejected on similar 
grounds. The V-K. color (1.87 ±0.04; iCutri et al.ll2003l; 
iDroege et al.l [2007) is consistent with an unreddened 
G8V (±1 subtype). At a distance of 39 pc, t he observed 
V magnitude (V = 10.50 : IDroege et~aLll2007f ) would cor- 
respond to My = 7.52, which is ^2 magnitudes fainter 
than typical G8V dwarfs. The colors of both candidates 
are also far too red to correspond to unre ddened young 
(~10 7 ~ 8 years) degenerate stars (see e.g. iKalirai et al.1 
12003! ) . Reddening is unlikely to be responsible for the 
brightness deficits of these two obje cts: reddeni ng is neg- 
ligible within 75 pc of the Sun (e.g. lLerov]|1993l), and HR 
8799 itself has negligible reddening (jZerbi et al.l 11999). 
We conclude that both of the two strongest NOMAD 
candidates identified bv IClose & Males (2009) arc prob- 
ably background field stars, unassociated with HR 8799. 

5.4. Search For Co-Moving Stars Within Ten Degree 

Radius 

Some of the famous young nearby stars in the solar 
vicinity have been found to belong to stellar associations 
which help in their age-dating (e.g. TW Hya, AB Dor), 
mostly through combining X-ray an d astrometry data to 
searc h for active, co-moving stars (Zuckcrm an fc Songl 
2004). This motivated us to search for evidence of a 
group in the vicinity of HR 8799 similar to other studies 
above. At present, we focus on X-ray /astrometric selec- 
tion since astrometric selection alone is rather daunting 

8 http:/ /www. projectpluto.com/photomet. html 
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and produces a subs tantial chance of false positives (e.g. 
iClose fc Maies|[2009l ). We searched a 10° (-10 pc) ra- 
dius around HR 8799 for signs of co-moving low-mass 
stars with a cross- r eferen ced list of X-ray-emitting stars 
(|Voges et all [1991 [2000l) in the Tych o -2 and UCAC2 
astrometric catalogs (|Hog et all 120001 : [Zacharia s et al.l 
[200l . As the 3D velocity of HR 8799 is well-defined, 
we tested whether the motions of these stars were con- 
sistent with motion towards HR 8799's convergent point 
(a, 6 = 99°.5 ± 2°.8, -31°. 1 ± 1°.8 deg). The most 
interesting system identified was the dMe star [K92] 94 
(= 2MASS J23081213+2126198J3 situated 1235" away 
from HR 8799. The star was o rigina lly discovered due to 
modest Ha emission bv iKunl (|1992|), and later classifie d 
as MOV without Ha emission bv lMartin fc Kunl (|1996[ ). 
However, if placed at 39 pc, its absolute K magnitude 
(M K ~ 8.1; based on 2MASS K s photometry) would be 
~3 magnitudes fainter than a typical young disk MOV 
star (M K ~ 4.8: lLeggettlll992f ) . Combined with the fact 
that its tangential motion differs from that of HR 8799 
by —4 km/s, it seems unlikely that [K92] 94 is associated 
with HR 8799. 

We therefore conclude that there still is no evidence 
for any co-moving stellar companions to HR 8799 which 
could help constrain its age. For the purpose of estimat- 
ing planet m a sses w e thus use the same age estimate as 
IMarois et al.l (f2008h of 60±*° Myr. 

6. COMPARISON TO GIANT PLANET MODELS AND 
BROWN DWARFS 

The observational results above can be used to con- 
strain the physical parameters for the planets around HR 
8799. We have used two separate approaches: compar- 
ison with young, giant planet models in the literature, 
and comparison with field brown dwarf photometry and 
spectra. 

For the giant planet models, we can use the evolu- 
tionary cooling tracks of lBurrows et al.l (|1997t ). hereafter 
B97, to estimate the size of the planets. The modeled 
radii of objects in the mass range of 5-20 Mj up and 
an age range of 30-160 Myr vary from 1.15-1.29 Rj up . 
We adopt a radius of 1.22 Kj up to estimate the flux den- 
sity from the model SEDs. This size suggests a gravity 
of log(g)=4.0-4.6 is expected for these objects. We use 
models of log(g)=4.5 for comparison with the observa- 
tions. 

Because of the small variations in radius versus age 
and mass, the key observable to constrain is the effective 
temperature of the objects. Depending on the age of 
the system, the effective temperatures will correspond to 
different mass estimates. 

6.1. Giant Planet Models 

Models of brown dwarf and giant planet atmospheres 
have been carried out with increasing levels of de- 
tail over the past decade. Non-g r ay atmosph eres have 
been comput e d by Burrows et al.l (|1997L 120031 ). and by 
iBaraffe et all (|2003[) that cover this temperature range. 
The effects of both clouds and non-equilib r ium c hem- 
istry have been studied by ISaumon et al.l (|2006f) and 



9 The originally published position from Kun ( 1992) is appar- 
ently off by ~30". We correctly matche d the 2 MAS S source to the 
star using the finder chart in Fig. 8 of Kun (1998) and SkyView 
( http: / / skyview.gsfc.nasa.gov/ ) . 



iHubenv fc Burrows! (|2007f ) . iHubenv fc Burrowsl (|2007l ) 
have developed a grid of models using various diffusion 
coefficients, reaction rates for CO — > CH4 and gravities. 
They find that the M band flux, in particular, is af- 
fected for T e ff ~ 1 000 K, exactly the range fo r the three 
HR 8799 planets. ISaumon fc Marlevl (|2008l and refer- 
ences therein) have developed models that incorporate 
vertical mixing (ISaumon et al.ll2006l; ILeggett et al.ll2007D 
and increased metallicitv (|Fortnev et al.ll2008D tomatch 
spectra of brown dwarfs and predict the SED's of giant 
planets. We compare the ob s ervati ons above to families 
of models bv ISaumon et all (|2006l ) (hereafter S06) and 
IHubenv fc Burrowsl (|2007l ) (hereafter HB07) that include 
these effects. 

We plot a color vs. magnitude at L' and M in Figure 
[2] to compare the observations with the models. The 
M band limits are the three a limits listed in Table 
03 The data are compared to both models that are in 
chemical equilibrium and models which incorporate non- 
equilibrium chemistry. While b could be fit by any of 
the model sets, the location of c is inconsistent with 
the equilibrium S06 models, and d is inconsistent with 
both sets of equilibrium models. For both the HB07 
models and the S06 models, the non-equilibrium se- 
quence uses an eddy diffusion coefficient (K 2Z ) of 10 4 
cm 2 s -1 . The HB07 models shown use the "slow" re- 
ac tion for the CO —¥ CH 4 c hemical reaction described 
in IHubenv fc Burrowsl (|2007[ ). The S06 models use an 
equivalent rate to the "fast2" rate described in HB07 
(D. Saumon, private communication). , The models in- 
corporate clouds and are calculated for log((7)=4.5 The 
temperatures for each sequence are in 100 K increments, 
from 800 K to 1600 K for the HB07 models and 800 to 
1300 for the S0 6 models. The color s suggest that, as with 
brown dwarfs (jLeggett et al.ll2007l ). the HR 8799 objects 
are best fit by models which incorporate vertical mixing 
into the atmospheres. 

Figure [3] shows a color-color plot for the object and 
models, using L'-M vs. [3.3]-L. The models shown are 
the same as for Figured] Due to its faintness, the limits 
on b are consistent with any of the models. The [3.3]-L' 
color of c and d is bluer than predicted by the equilibrium 
models over the range of temperatures considered. The 
HB07 non-equilibrium models are only consistent with 
the [3.3]-L' color for effective temperatures above 1500 
K, for c, and 1430 K, for d. The S06 non-equilibrium 
models are well matched to colors, for c, which corre- 
spond to 1200-1350 K and, for d, 1050-1250 K. For b, 
the effective temperature estimated from the S06 model 
colors is below 950 K. In general, the S06 models indi- 
cate cooler temperatures by ~300 K for all three planets, 
compared to the HB07 models. 

Model fits to the [3.3]-L' color can be used to estimate 
the effective temperatures of the planets. These fits are 
shown in Figure |4] The HB07 models result in T e g< 
1310 K for b, and T off =1620 and 1490 K for c and d. For 
the S06 models, the corresponding results are <940 K 
for b, and 1300, and 1170 K for c and d, respectively. 

The model fits to the colors do not constrain the overall 
flux. If a size of 1.2 Rj up is used to estimate the flux, 
consistent with expected size of these objects from the 
B97 models, the results are significantly brighter than 
the observed photometry at L' . To match the observed 
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flux at the best fit temperatures, the planets would need 
to be smaller than expected from the B97 models. For 
the HB07 models, planets c and d would need to have 
radii which are 44 and 63% of their expected size of 1.2 
Rj«p- For the S06 models the radii would need to be 56 
and 84% as big. 

Figure M again shows the SEDs of the HB07 and S06 
models compared to the HR 8799 photometry, but now 
matched to the absolute magnitude at L' (M L ). For each 
planet, a model temperature is selected that best fits the 
value of M^, using an assumed radius of 1.22 Rj up for 
the planets. The effective temperature and model pa- 
rameters corresponding to the fit for each plot is inset in 
the graph. For the non-equilibrium SEDs, the average 
modeled flux in each passband is plotted as a horizontal 
line spanning the passband. The difference in the best fit 
SEDs between the HB07 and S06 models suggests that 
the details of how the vertical mixing is modeled is im- 
portant for accurate estimates of colors, especially in the 
3.3 /im portion of the spectrum. 

Different diffusion rates, reaction rates and gravities 
,as well as a comparison of cloudy and clear models were 
explored for the HB07 models. These comparisons are 
shown in Figure |6] Only cloudy models are available for 
the various reaction rates tried, while the varying diffu- 
sion rates are only calculated for clear models. While 
this leads to an imperfect comparison, the general trend 
for each parameter can be seen. Cloudy models can re- 
produce the [3.3]-L', colors seen for c and d, but only for 
T e ff>1400 K. The other parameters do not significantly 
affect the [3.3]-L/ colors. 

6.2. Comparison to Brown Dwarf Spectra 

Another useful comparison to the observational results 
can be found in the colors of brown dwarfs in the field. 
Although field brown dwarfs of similar temperature to 
the HR 8799 are likely more massive, due to being older, 
they might still be expected to provide a good approxi- 
mation of the expected SEDs of these objects. 

Although the [3.3]-L' colors are blue relative to most 
model SEDs, similar colors are seen in spectra for brown 
dwarfs at these effective temperatures. In addition, the 
observed [3.3]-L' color changes significantly for spectral 
types in the L9 to T3 range, confirming that this region 
of the spectrum is useful for constraining the physical pa- 
rameters of brown dwarf and giant plane t atmospheres 
in thi s temperature range. For example, Lcggett et al. 

( 2008) analyze HN Peg B in detail and iStcphcns et al. 

( 2009) present infrared spectra of a large sample of brown 
dwarfs, including, 3-4 fim spectra of a subset of the ob- 
jects. The objects SDSS 1207+02, SDSS 0758+32, HN 
Peg B, and 2M 0559-14 with spectral types of TO, T2, T3, 
and T4.5 have [3.3]-L' colors of approximately 0.75, 1.2, 
1.9, and 2.5 respectively. HR 8799 c and d have a [3.3]-L/ 
color of 0.63 and 1.16 respectively. This suggests that the 
[3.3]-L' color of c and d are consistent with a T0+1 and 
T 2+1 dwarf, respec t ively. From the te mperature scale 
of iGolimowski et all (|2004| ) as revised by iStephens et al.l 
(2009) this would indicate T eff = 1200- 1300 K, for these 
objects, consistent with the S06 models, b has a [3.3]-L' 
color of > 1.9 constraining its spectral t ype to be later 
than T 2, or T cff < 1200 K. Of course, the [Stephens et al.l 
(2009) T c ff vs. spectral type relation shows scatter of 
~100 K about a given spectral type, suggesting caution 



should be used in comparing these results, 

Finally, we note th at the H-[4.5] c olor-te mperature re- 
lation developed by iWarren et al.l (|2007| ) can be used 
to constrain the temperature of these objects, assum- 
ing they have similar emerge nt spectra to field brown 
dwarfs. iLeggett et al.1 (|2007| ) derive a constant [4.5]- 
M'=-0.3 color co rrection for objects i n this temperature 
range. From the iMarois et al.l ([2008) H band photome- 
try we calculate an H-[4.5] color of less than 2.47, 1.63, 
and 1.56 for b, c, and d respectively. Using the Warren 
scale we derive temperature greater than 900, 1100, and 
1150 K for the objects, consistent with the color-derived 
estimates from the S06 models. 

7. DISCUSSION 

The key feature of the 3-5 (im results is the blue 
color of the photometry, compared to equilibrium mod- 
els. While the L'-M colors can be matched with the HB07 
models, the [3.3]-L' model colors are too red unless effec- 
tive temperatures of >1400 K are used. Vertical mixing 
can explain both the M band and 3.3 /im observations 
in the S06 models. The difference between these similar 
model families suggests that the details of how the non- 
equilibrium chemistry is modeled a ppears to be an im - 
portant factor in the resulting SED. Jan son et al.l ([2010D 
reach a similar conclusion from their analysis of c s spec- 
trum, in the L band. 

Even the S06 models have a discrepancy between the 
effective temperature derived from the [3.3]-L' color 
and the best fit according to the L' photometry, with 
the T e ff predicted to be ~ 100-400 K lower from the 
L'photometry compared to the color-derived tempera- 
ture. The correspondence of the color-derived tempera- 
tures from S06 models with those of field brown dwarfs 
suggest that the model temperatures are reasonable. 
Based on the effective te mperatures, the deriv ed masses 
from the cooling models of lBaraffe et al.l (2003) are 12+2 
Mj„ p for c and 11+2 Mj up for d, for an assumed age 
of 60 Myr. The temperature limit on b results in a mass 
limit of < 9 Mj up . The se mass estimates ar e consistent 
with the estimates from IMarois et ail (|2008l ). 

The large age range for the HR 8799 system (30-160 
Myr) increases the uncertainty in the estimated mass of 
these objects. This error term is larger than that from 
the photometry, resulting in a mass range of ~ ±3 Mj up 
for the objects. Combined with the photometric uncer- 
tainty, this indicates a mass of 12+4 Mj up for c, 11+4 
for d, and an upper limit on b of < 12 Mj up . 

The dimness of c and d, given the color-derived tem- 
peratures is problematic for any evolutionary model of a 
giant planet. Even a very old object will not be small 
enough to match the L' apparent magnitudes. The best 
fits to T G ff of c and d result in radii for the objects of 
0.7 Rj up and 1.0 respectively. Given the young age of 
HR8799, a possible explanation of this discrepancy is 
that dust extinction of approximately 1.3 and 0.4 mag- 
nitudes at L' is reducing the observed flux for c and d 
respectively. Since HR 8799 is relatively nearby an shows 
no sign of foreground extinction, this dust obscuration 
would need to be intrinsic to the objects themselves. 

Alternatively, lower gravity, metallicity, details of the 
chemical abundance versus height, or even incorrect 
opacity estimates for methane may be affecting the col- 
ors. Models with log(g)=4 for the S06 were also com- 
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pared to the observations, but were not different from the 
log(<?)=4.5 models shown. However, it is possible that 
gravity, combined with effects not properly captured in 
these models (such as metallicity or details of the chem- 
ical abundance variation versus height) may be affecting 
the colors for these objects in a more significant way than 
for field brown dwarfs. Typical field brown dwarfs have 
median ages of ~3 Gyr (jDahn e t al. 2002). Hence, the 
typical gravity of a T e ff= 1000 K object in the field is 
log(p) = 5.0, while these objects are expected to have 
log(g) = 4.3. If the vertical mixing is similar for both 
the field brown dwarfs and the planets, it is reasonable 
to expect that the stronger than expected CO absorp- 
tion, and the weaker than expected CH4 absorption is 
an effect of the lower gravity. 

Interestingly, the NIR colors of these objects, as well as 
2M1207b, are redder than field brown dwarfs of the same 
brightness (M arois et al.ll2008[ ) . The NIR color difference 
could well be explained by their youth, and thus their 
lower gravity. I ndeed, the young field L dwarf 2M0141- 
4633, found by (jKirkpatrick et al.1l2006|) . has nearly iden- 
tical NIR photometry to c and d,and is discrepant with 
other field L dwarfs. The estimated mass of 6-25 M j up for 
this object indicates it should have similar gravity to the 
HR 8799 planets. This suggests that the red NIR col- 
ors are a low gravity effect. Kirkpatrick et al. attribute 
the red colors to lower collision-induced absorption of 
H2, which primarily suppresses the K band. At lower 
gravities the H2 CIA is reduced, causing the redder H-K 
color. Thus, while the color trend is opposite in the ther- 
mal IR, it appears that both effects might be explained 
by the youth, and low gravity, of the objects. 

While in principle, the balance between CO and CH4 
in the se atmospheres is go verned by Le Chatelier's prin- 
ciple (Bu rrows et al.ll200lD . the details of how the verti- 
cal mixing occurs, appears to be important in predicting 
the emergent spectra. An addi tional effect may be the 
metallicity of the atmospheres ([Fortnev et al.lr2008f) . In 
higher metallicity objects, the formation of CO might 
expected to be enhanced since more oxygen is available 
relative to hydrogen, compared to lower metallicity at- 
mospheres. This suggests that giant planet atmospheres 
would be blue in the 3-5 /im region for more metal-rich 
objects. 

Altho ugh HR 8799 is me tal-poor star, detailed obser- 
vations (jGrav fc Kavell999l) indic ate solar abundance for 
lighter elements (<Sadakanell2006l) . a consequence of the 
A Boo phenomenon. This suggests that models with so- 
lar metallicity are probably appropriate for interpreting 
the colors of the HR 8799 planets. The solar metallicity 
(or perhaps metal-poor if the A Boo phenomenon is dis- 
counted) cannot explain the blue 3-5 /jm colors. If the 
formation of the planets caused them to be significantly 
metal-enhanced relative to their stars, this may be able 
to explain their colors. Further observations may be able 
to disentangle the metallicity of these objects from the 
other parameters which govern their SEDs. At the same 



time, detailed models are needed to explore the poten- 
tially complex interplay between gravity, vertical mixing, 
and metallicity. 

It is interesting to note that mass estimates, derived 
from model fits to the L' photometry are not significantly 
different from color-derived estimates. For example, the 
non-equilibrium fits from the HB07 models, for an age 
of 60 Myr, result in mass estimates of 9, 10.5, and 11 
Mj up for b, c, and d, respectively. If we use the S06 
models, the estimates are 7, 9, and 10 Mj up respectively. 
The HB07 are likely a better estimate of the effective 
temperature from the Movables, since these models treat 
the temperature profile in a self-consistent way with the 
vertical mixing, while the S06 models do not take this 
into account. The SED's for these objects, shown in Fig 
[SI indicate that the non-equilibrium are consistent with 
the L'-M colors, while being discrepant with the [3.3]-L' 
colors. Although the internal inconsistencies point to- 
ward a need for better models and further observations, 
the current results all indicate the objects are massive (> 
10 Mj up ), providing a challe nge to models for explain- 
ing the stability of the system (iFabrvckv fc Murrav-Clavl 
I2008t [Reidemeister et al.ll2009f) . as well as its formation. 

8. SUMMARY 

Measurements of HR 8799 in thermal IR have con- 
strained the SED of these planets, helping to better 
understand their physical parameters. The L'-M and 
[3.3]-L'colors of the objects appears to be blue, relative 
to equilibrium models. Models which account for non- 
equilibrium chemistry of CO and CH4 can account for 
this effect, but the details of the modeling leads to dif- 
fering results, depending on the input assumptions. Im- 
proved modeling of these objects will help understand 
the source of this. Deeper 3.3 ^m and M band photom- 
etry will be helpful in constraining the SEDs of these 
interesting objects. 

The effective temperatures of 6, c, and d, derived from 
the colors as predicted by the S06 models are <940, 1300, 
and 1170 K, indicating the objects are < 9 Mj up for b and 
12, and 11M j U p for c and d for an assumed age of 60 Myr. 
Th ese estimates provide an independent constraint from 
the IMarois et al.l (|2008f ) results on the effective temper- 
atures and masses of the companions to HR 8799. The 
results confirm the massive nature of this planetary sys- 
tem, and highlight the need for more detailed models to 
fit the colors of these objects. 
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Fig. 1. — Images at 3.3/xm, L' (3.8 fim) and M (4.8 fim) of HR 8799. All three components are detected at L' (center). At 3.3 ^tm(left) 
only c and d are detectable. At M band (right) none of the planets are detectable above the sky background noise. The circles identifying 
the locations of the planets in each image are the approximate size of the photometric aperture used. 
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L'-M 

Fi g. 2. — Color magnitu de relation for models compared to observed photometry. Model colors are shown from Hubcny & Burrows (2007) 
and lSaumon et all H2006T ) . The temperature range spans 800-1600 K for the HB07 models and 800-1300 K for the S06 models. The models 
use K zz =10 4 cm 2 s — 1 , log(g)=4.5, and a cloudy atmosphere model. A chemical reaction rate corresponding to the "slow" rate is used for 
the HB07 models, and the "fast2" rate is used for the S06 models. 



12 



Hinz et al. 




Fig. 3. — Color-Color plot for L'-M vs. [3.3]-L'. Models from lHubenv & Burrows! H2007T ) and lSaumon et al.1 l|2006l ) are shown. The same 
parameters are used as in Figure 2. 
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Fig. 4. — Comparison of HR 8799 photometry to models by[Hubcnv & Burrow 1 (200^ and ISaumon et alJ l(2006h for 6(top), c(middle), 
and d(bottom). For each planet, models are shown that best match the [3.3]-L' color. The overall flux is allowed to vary to fit both the 
[3.3] and L'photometric points. The band-integrated flux density of the models are shown as horizontal lines for the observed passbands. 
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Fig. 5. — Comparison of HR 8799 photometry to models by Hubcnv & Burrows (2007) and Saumon et al. (2006) for 6(top), c(middle), 
and d(bottom). For each planet, models are shown that best match the L' band photometry. The chemical equilibrium models appear to 
over-predict the M band flux, while non-equilibrium models are consistent with the 3 a upper limits. The band-integrated flux density of 
the non-equilibrium models are shown as horizontal lines for the observed passbands. 
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Fig. 6. — Color-Color plots using the Hubcny & Burrows (2007) models. The upper left compares different diffusion rates. The right 
compares different reaction rates. The lower left compares varying gravity, while the lower right compares cloudy and clear models. The 
models are unable to reproduce the [3.3]-L'colors unless temperatues above T e ff=1400 K are used. 
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TABLE 1 

Log of data acquisition for Clio observations of HR 8799 

Date Total Integration (s) Filter Total Rotation (degrees) 

21 November 2008 5694 L' 72.0 

21 November 2008 9600 Barr M 31.8 

9 January 2009 2907 3.3 fim 5.3 

12 September 2009 6780 3.3 (im 30.2 



TABLE 2 
Adopted Photometric Values for 
HR8799 



Band 


Stellar Magnitude 


zeropoint (Jy) 


3.3 fim 


5.23 


330 


L' 


5.22 


235 


M 


5.21 


154 



TABLE 3 

Photometric Results 



System 


Band 


b 


c 


d 


Apparent Vega Magnitude 


3.3Atm a 
L' b 

M c 


>1 7.8 
15.81±0.23 
>14.7 


15.63±0.3 
14.99±0.12 
>14.7 


15.53±0.3 
14.39±0.19 
>14.7 


Absolute Vega Magnitude 


M3.3 
M L , 
M A / 


>14.82 
12.83±0.23 
>11.72 


12.65±0.3 
12.01±0.12 
>11.72 


12.55±0.3 
11.41±0.3 
>11.72 


Flux at 10 pc 


F3.3 (mJy) 
F L , (mJy) 
F M (mJy) 


<0.39 
1.73±0.41 
<3.3 


2.87±0.86 
3.69±0.44 
<3.3 


3.15±.95 
6.41±1.9 
<3.3 



a Custom filter with zeropoint listed in Table [2] 

MKO photometric system ( Tokunajja & Vacca 2005). 
c Bcsscll and Brett photometric system (Bcsscll & Brett 1988). 



TABLE 4 

ASTROMETRIC RESULTS FOR COMPANIONS TO HR 8799 



Date 


Band 


HR 8799 b 


HR 8799 c 


HR 8799 d 






E ["], N ["] a 


E["],N["] a 


E["],N["] a 


21 November 2008 


V 


1.542±0.01, 0.780±0.014 


-0.631±0.015, 0.671"±0.02 


-0.215±0.021, -0.644±0.013 


9 January 2009 


3.3 fim 




-0.612±Q.03, 0.665" ±0.03 




12 September 2009 


3.3 fim 




-0.625±0.02, 0.725" ±0.02 


-0.282±0.03, -0.590±0.03 



a Offset from HR8799. 



